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In supersymmetric models of warm inflation, the large temperature of the radiation bath pro- 
duced by the dissipative motion of the inflaton field may induce a significant thermal abundance 
of potentially dangerous gravitinos. While previous discussions of this problem focused on grav- 
itino production only at the end of warm inflation, similarly to conventional reheating scenarios, 
we study the full evolution of the gravitino abundance during and after inflation for simple mono- 
mial potentials, taking into account the enhanced gravitino and potentially gaugino masses due 
to supersymmetry breaking during inflation and the smooth transition into a radiation-dominated 
era. We find, on one hand, that the continuous thermal production increases the gravitino yield, 
although, on the other hand, 'freeze-out' occurs at temperatures much lower than previously esti- 

^N) mated. Moreover, for sufficiently strong dissipation, which allows for sub-planckian inflaton values, 

( the lower radiation temperature significantly alleviates and potentially solves the gravitino problem, 

with a baryon asymmetry being nevertheless produced through dissipative effects. Our analysis may 

04 also be relevant to standard reheating as an oscillating inflaton will also change the gravitino mass, 

potentially modifying the produced gravitino yield. 
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I. INTRODUCTION 



< 

I Inflation [B-S] has been incredibly successful in providing solutions to the problems of the standard cosmological 

model. It can set the initial conditions, which give rise to the high degree of flatness and homogeneity that we 
^ observe in the universe today. From particle physics motivated models, it not only yields a mechanism for accelerated 
expansion but also explains, through quantum fluctuations, the origin of the temperature anisotropies in the Cosmic 
Microwave Background and the seeds for the observed Large Scale Structure. 

In the standard cold or isentropic inflation scenario, the early universe is dominated by the vacuum energy of a 
scalar field which is slowly rolling down its potential, resulting in a period of accelerated expansion. This occurs whilst 
its kinetic energy is negligible compared to the potential energy until, at some point, the potential typically steepens 
and the inflaton begins oscillating about the minimum of its potential. During the period of accelerated expansion, 
the inflaton is assumed to have negligible couplings to other fields in order to keep the potential flat enough for a 
sufficient number of e-folds of inflation (~40-60) to occur and as a result the universe supercools. However, once it 
begins oscillating, there must be interactions that convert the vacuum energy into radiation in order to reheat the 
universe. 

Hence, the inflaton cannot be an isolated system and, while the standard picture assumes that any interactions 
have a negligible effect on the dynamics of inflation and only become important during reheating, this need not be 
the case. In the alternative warm inflation paradigm, such interactions may in fact lead to dissipation of the inflaton's 
kinetic energy into light degrees of freedom, which in the simplest case may thermalise, resulting in the presence of a 
^ nearly-thermal bath concurrent with the accelerated expansion. 

In the early universe, gravitinos can be abundantly produced, potentially leading to overclosure of the universe or 
spoiling the abundances of light elements predicted by the standard big bang nucleosynthesis model (BBN). While 
in cold inflation thermal production of gravitinos occurs only during the reheating phase, in warm inflation this 
is concurrent with inflation due to the presence of a thermal bath. In standard reheating, gravitino overproduction 
constrains the reheat temperature, i.e. the maximum temperature after inflation when the universe becomes radiation 
dominated. There is, however, a certain amount of tension in this case between having a large enough reheat 
temperature to allow for a thermal mechanism for baryogenesis, whilst keeping it low enough to avoid overproducing 
gravitinos. In warm inflation, on the other hand, this tension can be relieved, as a baryon asymmetry may in fact be 
produced at low temperatures through dissipative effects |4], potentially avoiding overproduction of gravitinos. 
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Gravitino production in warm inflation has been considered previously in [51 [B] , where it was assumed that the 
effective 'reheat temperature' occurs when the radiation energy density becomes equal to the inflaton energy density, 
P<ji — PRi E^J^d that standard reheating constraints on gravitino production can be applied. This may, however, 
overestimate the temperature at which the gravitino yield freezes out, as radiation does not yet fully dominate the 
energy density at this stage. Moreover, standard reheating constraints may not a priori be applied in warm inflation 
scenarios due to the non-negligible abundance of gravitinos produced during inflation, which may potentially lead to 
a larger yield. This may, in fact, be the case also in conventional models, since the reheating phase is not necessarily 
instantaneous and thermal production of gravitinos may potentially occur for the duration of reheating and not freeze 
out until the universe is fully radiation dominated, resulting in a cumulative effect similar to that of warm inflation. 
Finally, we note that supersymmetry is broken during inflation, leading to gravitino masses parametrically close to 
the Hubble parameter and potentially to massive gauginos, which may also modify the production rate during warm 
inflation. Similarly, this may change the standard reheating constraints, as the gravitino mass also varies during the 
oscillating phase. 

With these new insights in mind, we revisit the production of gravitinos in supersymmetric warm inflation, numer- 
ically evolving the Boltzmann equation for gravitinos into the radiation era. In Section [ll] we give a brief review of 
the standard gravitino cosmology in cold inflation and in Section III we outline the basic features of warm inflation, 
focusing as a working example on monomial potentials in the sub-planckian regime. We discuss thermal gravitino 
production in warm inflation in Section |IV| and present results for stable and unstable gravitinos, considering the 
effects of inflaton-dependent gaugino masses in both cases. In Section [V] we summarize our main results and discuss 
possible directions of future research in this topic. 



II. STANDARD GRAVITINO COSMOLOGY 



Supersymmetry is an attractive theory for inflationary dynamics due to the presence of a whole host of scalar fields, 
as for example the superpartners of Standard Model quarks and leptons, automatically protecting the scalar potential 
from quadratic loop corrections that may spoil its required flatness. However, several single field models typically 
require inflaton expectation values close to the Planck scale, where supergravity effects start playing an important 
role. The gauge particle of supergravity is the massless spin-3/2 gravitino and, when supersymmetry is broken, 
the gravitino becomes massive and absorbs the spin-1/2 goldstino through the super-Higgs mechanism. Due to its 
indiscriminate coupling, the neutral gravitino couples to all fields universally, whether in Standard Model/visible sector 
or other hidden/sequestered sectors, with planck-suppressed interactions, making it a potential candidate for dark 
matter. This suppressed coupling makes it, however, unlikely to be detected at man-made colliders. Gravitinos may 
nevertheless be abundantly produced in the early universe through a variety of thermal and non-thermal processes 
due to the large energies involved. Unfortunately, due to our ignorance of the mechanism behind supersymmetry 
breaking, its mass is unknown and can only be constrained by cosmological considerations. 

Without a period of inflation the constraints on the gravitino mass are quite severe. In the standard cosmological 
model, the early universe is radiation-dominated and, at early enough times, the temperature will be high enough for 
gravitinos to be in thermal equilibrium with the radiation bath. For stable gravitinos, the early freeze-out associated 
with planck-suppressed interactions can thus result in overclosure of the universe unless m^, ^ 1 keV . For an 
unstable gravitino, its mass needs to be larger than ^ 10 TeV, otherwise it will decay during BBN and spoil the 
predictions for the light element abundances [H]. With a period of inflation, any initial population of gravitinos is 
diluted away and, in cold inflation models, no gravitinos are thermally produced until reheating. These strict bounds 
on the gravitino mass are thus somewhat relaxed and replaced by upper limits on the reheat temperature. 

Gravitinos are primarily produced by the scattering of particles in a thermal bath. Due to the stronger coupling, 
the dominant production comes from inelastic 2 — )■ 2 QCD processes involving left handed quarks (5), squarks (g), 
gluons (5), gluinos {g) and gravitinos (G) such asg + g^gA-G, q + q g + G and q + g ^ q + G. While this 
contribution to the thermal gravitino production rate in supersymmetric QCD has been calculated in [5], in this work 
we will adopt the complete SU{3)c x SU{2)l x U(1)y thermal production rate computed in [10]. Gravitinos can also 
be produced from the decay of the inflaton during its oscillating phase after inflation (see [TT| - fT7] V However, as we 
discuss below, this only becomes significant during the radiation era and, due to the suppression of inflaton oscillations 
for strong dissipation, it is subdominant for the monomial models we consider. They can also be produced through 
the decay of other particles in the thermal bath, assuming this is kinematically allowed. However, during inflation 
the Hubble parameter is in general much larger than the relevant planck-suppressed decay widths, H ^ r^ecai/i so 
that these decays will also not occur until the Hubble parameter drops signiflcantly in the radiation era. As it is our 
aim to highlight the differences between the gravitino production in warm and cold inflation, we will thus focus on 
thermal production processes. 

If gravitinos are abundantly produced in the early universe, this can pose a problem for inflationary model building. 
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There are two situations to consider in gravitino cosmology, depending on whether the gravitino is the Hghtest 
supersymmetric partner (LSP) and stable by virtue of R-Parity conservation, or otherwise the gravitino is unstable 
and will decay at some stage in the cosmological evolution. Constraints on primordial abundances come in either case 
from two sources. Firstly, the abundance of the LSP must not exceed the observed dark matter abundance ^ [TS]: 

r!z5M/i2 - 0.105™ . (1) 

Secondly, the decay products of the next to LSP (NLSP) must not spoil BBN predictions for light element abundances. 
Radiative decay of the NLSP where photons and charged particles are emitted can induce electromagnetic showers, 
disintegrating the light elements. The NLSP can also decay into quarks or giuons, which subsequently hadronise. These 
hadrons can then induce interconversions between the background protons and neutrons, enhancing the neutron to 
proton ratio and thus resulting in an overproduction of '^He. The energetic nucleons can also destroy the background 
"'He and non-thermally produce D, T, ■^He, ^Li and ^Be. There is also the possibility that if the LSP is charged then 
it could bind with background nuclei and change the nuclear reaction rates, in particlar that of ^Li. 

If the gravitino is unstable and m^, ^ 20TeV, its lifetime is longer than 1 s [12] and so it will be subject to the 
BBN constraints mentioned above. It will also decay into the LSP which needs to satisfy the dark matter constraint 
in Eq. ([T]) . If the gravitino is stable then it must satisfy the dark matter constraint ([T]) and the decay of the NLSP 
into the gravitino must avoid upsetting BBN predictions. For more details on BBN constraints on LSP and non- 
LSP gravitino primordial abundances see [I9H26] . The problems mentioned above constitute the so called "gravitino 
problem" . 

The number density of gravitinos is described by the Boltzmann equation: 

Tig + 3Hn^ = Ca , (2) 

where we neglect gravitino decay and gravitinos produced from decays of other fields. The collision term, Cg, describes 
gravitino production in a thermal bath and is given by |10j : 

The index i runs over the gauge groups [U [1)y , SU [2) l , SU {?>) c) where rug^ are the gaugino masses, gi are the 
gauge couplings and q = (11, 27, 72), ki — (1.266, 1.312, 1.271). The reheating phase is assumed to be instantaneous, 
immediately entering the radition era. Defining the gravitino-to-photon yield, Yq — n^/n-y, with = 2^(3)T^/7r^, 
and assuming that TR = constant, where R is the scale factor, we obtain from Eq. ([2]): 

^ = ^ ^4^ 

dT H{T)Tn^{T) ' ^ ' 

As Cq ~ r^, with only a mild temperature dependence from the couplings and gaugino masses, Eq. Q can be 
approximately integrated. Assuming that any initial population of gravitinos before the reheating phase is diluted 
away, Yq[Tr) — and that we are interested in the yield of gravitinos at temperatures T <^ Tr, e.g. at BBN, then: 

" HiTn)!^fiTn) ' 

However, TR only remains constant away from particle mass thresholds, as it is instead the entropy density that is 
conserved, sR^ ^constant. We can take this into account by diluting the yield: 

yoiT.) = 'W/rW^'^oiT.) = ^r^(T.) , (6) 
s{T2)/n^{T2) 9*[T2) 

where (T) is the number of relativistic degrees of freedom in thermal equilibrium at temperature T. Note that this 
is generically well below the gravitino yield in thermal equilibrium. We then obtain the following expression for the 
present abundance of gravitinos: 




, f W (7) 

gij \100GeVJ {lO^'^GeV J ^ ' 
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The subscript '0' indicates the present day value, with Tq = 2.73 K, .g*(To) = 3.91. We use the MSSM value for 
g*{TR) = 228.75 and the critical density Pc = 8.1 x IQ-^^/i^ GeV* with the constants = (0.018,0.044,0.117). It is 
understood that the couplings and masses should be evolved with the temperature. This provides the standard cold 
inflation constraints on the reheat temperature and gravitino mass to avoid overclosure for an LSP gravitino. It is 
evident that for m^, « 100 GeV, to avoid overclosure, the reheat temperature Tj^ ^ 10^° GeV. If the gravitino is the 
NLSP, then each gravitino will decay into one LSP and the primordial gravitino yield in Eq. ([s]) can be converted into 
the LSP yield through: 

nLSph' = "^np^ . (8) 

We can see that VL^^h^ > 0.105 is allowed as long as niLSP I'm-Q sufficiently small. For unstable gravitinos with 
rriQ ^ 20 TeV the strongest constraints come from BBN abundances, whilst above this the dark matter constraint 
for the LSP dominates [19.;. These constraints are often defined in the literature in terms of the gravitino-to-entropy 
yield, = Uq/s, which can be easily related to the more convenient definition in terms of the photon energy density 
used in Eq. We wih take the conservative bounds of Y? < IQ-^^ for 100GeV< < 1 TeV and Y? < IQ-^'^, 
^ 10~^^ for lTeV< < 3 TeV for branching ratios into hadrons of = lO"'^, B/i — 1 respectively [33]. For 
stable gravitinos the BBN constraints on the primordial yield from NLSP decays are quite model dependent, varying 
upon which particle is the NLSP, as well as its thermal abundance and mass. For more details on scenarios with 
sneutrino, slepton and neutralino NLSPs, see e.g. [20} [2T|. 



III. WARM INFLATION 



Warm inflation [37H35] (see also [3Uti32| ) is an alternative picture of inflation, where the inflaton has non- negligible 
interactions with other fields that lead to fiuctuation-dissipation dynamics and associated particle production con- 
current with accelerated expansion. While several efforts in the literature to analyze non-equilibrium dynamics fall 
within this category (see e.g. [33H36) ). the most well understood and extensively studied scenarios in the context of 
quantum field theory consider the case where radiation is produced in a nearly-thermalized state [291 l37H41j . This 
allows one to make accurate predictions for the effects of dissipation and has several attractive features from the 
model-building point of view. Firstly, the dissipative dynamics acts as an additional source of damping that allows 
for longer periods of accelerated expansion, which is particularly important in supergravity/string theories, where 
F-term supersymmetry breaking typically induces large inflaton masses, thus alleviating the associated eta-problem 
|42H46j . Secondly, whilst one assumes radiation to be a subdominant component of the energy balance in the universe 
for accelerated expansion to occur, in several scenarios it may actually come to dominate at a later stage, providing 
a smooth transition into a radiation-dominated era, in alternative to the standard reheating picture. Finally, since 
thermal fluctuations overcome the quantum vacuum fluctuations for temperatures T > the spectrum of primordial 
density fluctuations may be significantly modified [28j ,39j l47l - l49] , in particular suppressing the amplitude of tensor 
perturbations and inducing potentially observable deviations from a gaussian spectrum [50H53] . 

The dissipative dynamics arising from interactions of the inflaton with other fields arises through time non-local 
contributions to its quantum effective action, which for a slow-rolling inflaton may in general be computed using linear 
response theory. In the adiabatic regime where <j)/(j) < r"^, where t is the typical relaxation time of the nearly-thermal 
ensemble, this yields an effective friction term Tcf) in the equations of motion, which can be written as: 

4> + 3H{l + Q)(}) + V^ = , (9) 

where Q = T/3H and denotes the derivative of the potential with respect to the inflaton field. Noting that the 
effective density and pressure of the inflaton condensate are given by = 1^(0) andp^ = 4>'^ — V{(j>), respectively, 
this can be rewritten as: 

+ 3H{p^ + p^) = (^0) 

The energy lost by the inflaton field through dissipative effects is then gained by the produced particles (see e.g. 
[54j ) ■ and for g^, relativistic degrees of freedom this yields the following evolution equation for the radiation density, 
PR = ^2.g*TV30: 

PR + mpR = T4>\ (11) 

with inflation occuring for p^. This nevertheless allows for T > H, as mentioned above, in which case one may 

also neglect the quasi-de Sitter expansion when computing the dissipation coefficient in different quantum field theory 
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models. On the other hand, for T < H we expect the inflationary dynamics to be similar to the more conventional 
cold scenarios. 

Accelerated expansion occurs in the slow-roll regime, where V{(f>) ^ <j)^, (j) H(j). In warm inflation, this can be 
translated into the modified slow-roll conditions: 

e0 = ^(^)'<l + Q, m = ml(yp^^<l + Q, a, = ml^<l + Q. (12) 

In addition, we also require three more conditions. Firstly, we need the variation of T with respect to to be 
sufficiently slow, in order to avoid dissipation increasing too quickly and radiation dominating too soon. Secondly, 
we require that radiation is produced faster than it is diluted by the expansion of the Universe and, finally, that both 
quantum and thermal corrections to the infiaton potential are not too large and, in particular, do not induce a large 
inflaton mass. Once these conditions are violated either the radiation energy density starts to dominate or the inflaton 
is no longer overdamped and slow-roll ends (see e.g. |55j). In the slow- roll regime the equations of motion reduce to: 

Earlier attempts to construct models of warm inflation considered a direct coupling between the inflaton and the 
light fields that form the radiation bath, but in this case a sufficiently large dissipation coefficient also induces a large 
thermal mass to the inflaton field, which makes it difficult to achieve a sufficiently long period of accelerated expansion 
[371 156j . A more promising avenue considers a two-stage mechanism [40| , where the inflaton is coupled to heavy fields 
that may subsequently decay into light degrees of freedom. This is also a more natural approach since couplings 
to the inflaton generically induce large masses. Moreover, in supersymmetric models the leading corrections to the 
inflaton potential are logarithmic in this regime |57| , despite supersymmetry being broken by the finite temperature 
and energy density, keeping the flatness of the potential stable against quantum and thermal corrections. A generic 
superpotential implementing this mechanism is given by |551 158j : 

W = W{<^) + g^X'^ + hXY'^ . (14) 

The scalar component of $ is the inflaton field, with expectation value = ^p/\/2, which we assume to be real. Both 
the bosonic and fermionic components of the superfield X then acquire masses proportional to ip and can decay into 
the Y scalars and fermions, which remain light and form the radiation bath. For T <^ mx and a broad range of 
couplings and field multiplicities, the leading contribution to the time non-local effective action corresponds to 1-loop 
diagrams involving virtual X-scalars, and has been discussed in |58H60j . yielding a dissipation coefficient of the form: 



T«C^^, (15) 
where Ca, depends on the coupling h and the field multiplicities in the X and Y sectors. We restrict our analysis to a 



dissipative coefficient of the form (15). Although our analysis depends upon the form of the dissipative coefficient we 
expect our methodology to be applicable to other forms [371 [Ml 110] and our qualitative results to be similar. In this 
work, we will take as a free parameter of the model, bearing in mind that large values for this constant may require 
a somewhat large field multiplicity, which may be attained, for example, in GUT models with large representations 
or the multiple D-brane constructions described in [46| . 



A. Monomial potentials 



As a working example, we will take the inflaton superpotential to be of the form: 

Timr^ ■ 



(16) 



For r = 0, A < we recover supersymmetric hybrid inflation, with the X scalars corresponding to the waterfall 
fleld(s), and for r > 1 we recover chaotic inflation models. Assuming a canonical Kahler potential for the inflaton 
fleld, K(^, $^) — <I>^$, this results in the following scalar potential: 



V = X' 



M 



2r 



1 



1 

r + 1 



M 

nip 



2- 



r + 1 



1 

{r + iy 



M 

nip 



cxp 



(17) 
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Given our ignorance of fundamental quantum gravity effects, we will restrict our analysis to the sub-planckian 
regime |0| ^ rrip, where supergravity effects may also be ignored and the potential takes the simpler monomial form: 



The slow-roll parameters are given by: 



7]^ = 2r(2r- 1) 



£0 



2r 



-11 



2r - 1 



2r - 1 

From the slow-roll equations of motion (|13|), we can derive the following relation between Q and 

1/3 / \ / fj\^/^ 



where Cu = g*7r^/30. We can invert this to get: 



9C| Q{l + Qf 



(18) 



(19) 



(20) 



(21) 



where s = 1/(14 — 2r). The evolution of Q during inflation is found by differentiating Eq. (20 1 with respect to the 
number of e-folds. 



dQ 



Q 



(lOe^ - 6770 -I- 8a^) 



1 + 7Q 

Qi+^-(l + Q)i^- 
1 + 7Q 



2r / 9Cl 



(22) 



It is clear from Eq. (|22|) that for < r < 7, Q increases during inflation. The number of c-folds of inflation can then 

1 + 7Q 



be obtained by integrating Eq. ( 22 ) , giving 

'"3= dN, 



* dQ 



dQ = Cq 



Q, 



Qi+2«(i + g)i 



:dQ , 



with 



Co = 



2r 



(23) 



(24) 



The 'e' subscript denotes the number of e-folds at which the slow-roll conditions are violated, while the subscript 
indicates the value when cosmological scales leave the horizon during inflation. Performing the integral, this yields: 



where 



N, = CQ[Fr{Qe)-Fr{Q,)] , 



14a; 1 
^2i^i(l - 2a:, 12a;, 2 - 2a;, ~x) - -2Fi{-2x, 12x, 1 - 2a;, -a;) 



(25) 



(26) 



and 2Fi{a^ b, c, z) is the hypergeometric function. As Q — > oo, the number of e-folds approaches a constant and so it is 
not always possible to achieve the desired number of e-folds of inflation in areas of parameter space where Q diverges 
too early, corresponding to the breakdown of the slow-roll approximation. 

It is important to ensure that T > if, so that the dissipative coefficient in Eq. (151 can be calculated neglecting 
expansion effects. If we set T»/if* > 1 at horizon crossing, then this will hold for the duration of inflation, for 
< r < 7, as we can see from Eqs. (21 1 and (27): 



T 
H 



9Q 

C^X^ 



1/3 



2(l-r)/3 



(27) 
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As mentioned earlier, choosing as our other free parameter, the slow-roll dynamics are fully determined. We can 
use the amplitude of the primordial power spectrum, R 



1/2 



n 



2^ 



ml 



5 X 10-5 [ST, to fix 



(28) 



Eq. (20), we arrive at: 



Combining Eq. (28 1 with the form of the dissipation coefRcient in Eq. (15) and the relation between Q and </> in 

(29) 



1 2 



Note that only depends on {T^/H^) and not on or the form of the potential. Once we have Q*, we can integrate 
Eq. (23) to obtain the total number of e-folds. The regime where |0| <C mp corresponds to the strong dissipation 
limit, ^ 1. Making this approximation, we have: 
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1 



2rs VQ*" 



1 

yirs 



Cr 



(30) 



In Figure [T] we show the region of parameter space for the quartic (r = 2) and quadratic (r = 1) potentials where 
we can ignore supergravity corrections with a reasonable number of e-folds of inflation. We note that the quadratic 
potential, being flatter, requires lower values of than the quartic model to achieve the same number of e-folds of 
inflation in the sub-planckian regime. Notice that, although we need somewhat large values of to obtain 40-60 e- 
folds of inflation, this allows inflation to occur at sub-planckian field values, which is not possible in standard inflation 
and is therefore a very attractive feature of warm inflation. 




FIG. 1: Total number of e-folds for the quartic (r — 2, left) and quadratic (r = 1, right) potentials. The region where |(^| < mp 
is to the right of the labelled line. The shaded region corresponds to between 40 and 60 e-folds of sub-planckian inflation. 



For monomial potentials, we can also derive the following relation: 



1 + Q 



2(2r- 1) 



Q 



Q 



(31) 



We can thus see that, when the radiation energy density becomes equal to the inflaton energy density, PR — P(t,~ V , 
the slow-roll condition 77^ < \ + Q has already been violated. This means, in particular, the breakdown of the slow-roll 
equation for the radiation energy density, in Eq. (13), as pR becomes significant and radiation soon takes over. We 



wish to ultimately calculate the gravitino yield after inflation and this means that we need to evolve the full set of 
equations into the radiation era. To do this we must numerically solve the equations of motion ([9| and (11), which 
we will discuss in the next section. 
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IV. GRAVITINO PRODUCTION IN WARM INFLATION 



A. Particle Masses 



In the presence of supersymmetry breaking, the gravitino gains a mass: 

ttiq = mpexp(— G/2) . 



(32) 



For the monomial superpotential in Eq. (161 and a canonical Kahler potential, the gravitino mass is then given by: 

r+l 



Xnir, 



exp 



Ml 

2m2 



?' + 1 \ rrip ^ 

Comparing this to the Hubble parameter during inflation w {V/irrip), we get 

V3 f\cj)\ 



(33) 



niQ 
H 



1 V mr 



1 



1 



101 



r + l \rn 



2- 



r + 1 



1 



{r + iy 



-1-1/2 



As discussed above, we are interested in the sub-planckian regime, for which: 

_ 73 / ]0[ 
H r + 1 \ rriri 



(34) 



(35) 



In Figure [2] we can see that even if inflation is sub-planckian the gravitino mass can be a non-negligible fraction of 
the Hubble parameter during inflation, resulting in the gravitino mass being typically well above the TeV scale, in 
contrast with what was assumed in earlier works [6 . For example, with a quartic potential, where A ~ 10^^ yields 
the observed amplitude of density perturbations, if (p/mp Ri 1/2 then m^, ^ lO^*' GcV. 




FIG. 2: rriQ/H for various monomial potentials. 



As discussed earlier, thermal production of gravitinos procedes through the scattering of gauge bosons, gauginos, 
quark and squarks. It is, in particular, strongly dependent on the ratio of gaugino to gravitino masses, mg/rriQ. Having 
seen that supersymmetry breaking during inflation results in massive gravitinos, it is interesting to also consider its 
effect on the gaugino masses, which are given by the quadratic term in the Lagrangian: 

^gaugino = ^ exp(-G/2)G'(G-i)f ^A"A^ , (36) 

where fa/3 is the gauge kinetic function, which is a holomorphic function of the chiral superfields in the model. It is 
dimensionless and symmetric with respect to its two adjoint indices and, in renormalisable theories, it is proportional to 
Sap I g\ ■ Whether this function depends or not on the inflaton field is a model-dependent question and for completeness 
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we will consider both cases separately. Interesting examples of inflaton-dependent gauge kinetic functions may arise 
in extra-dimensional theories such as superstring/M-theory, where the inflaton is identified with a modulus field (see 
e.g. Considering the case where the infiaton's supersymmetry breaking effect is communicated to the visible 

sector through gravitational interactions, we can expand the gauge kinetic function in powers of (p/nip, yielding: 

/a/3 ~ f ^ + /q— + ... ) , (37) 
\9l J 

where fa is a dimensionless coupling, which for simplicity we will assume is universal to all the gauginos and will take 
to be 0(1). Although the infiaton field modifies the gauge couplings, this will not change the running of the couplings 
significantly since are considering sub-planckian field values. With the above expansion for the gauge kinetic function, 
the gaugino mass is given by: 

"~ M 

2m 




/-P(^)> (38) 



and for sub-planckian field values this reduces to: 



A \ I f^^f^^^f- (39) 
4 \ nip / Avip 4 

The gaugino masses are thus proportional to the Hubble parameter. We then find that: 

-^^^i±3i(±y\ (40) 

niQ 4 \ rup J 

so that gauginos are gencrically heavier than gravitinos during inflation, which will have important consequences on 
gravitino production. 

If the only source of supersymmetry breaking were the inflaton superpotential, then it is evident that as the infiaton 
rolls to its minimum supersymmetry would be restored. This is obviously not the case in nature, and so we will consider 
a supersymmetry breaking contribution from a hidden sector that gives rise to TeV-scale supersymmetric partners. 
The details of this hidden sector will not be important to the thermal production mechanism and so we can take the 
following phenomenological approximation for the masses: 



(41) 
(,2) 



where the subscript denotes the inflaton contribution and '0' indicates the low-energy hidden sector contribution 



B. Gravitino yield evolution 

We numerically solve the warm inflation equations Q and along with the Boltzmann equation for the gravitino 
number density, Eq. ([2]), with the collision term given by Eq. ^. We run the couplings and gaugino masses with 
temperature at one-loop assuming they unify at the GUT scale, Tqut = 2 x 10^^ GeV, with universal gaugino mass 
Tni/2 = 400 GeV. For convenience, we evolve the equations in terms of number of e- folds, Hdt — dNg, which we will 
use both during and after inflation. We find that the thermally produced gravitino yield freezes out and approaches 
a constant value after inflation ends when the following three conditions are met: 

• The gravitino has settled to its low-energy mass, given by the hidden sector contribution tHq^^ ; 

• The universe is in the radiation-dominated regime, where pR ^ exp(— 47Ve), i.e. radiation must have ceased 
being signiflcantly produced by dissipation; 

• Thermal production of gravitinos must have stopped, so that the collision term in the Boltzmann equation is 
negligible and thus the number density evolves as ~ exp(— 3iVe). 
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A helpful consequence of being in the sub-planckian regime is that the large value of makes the primordial yield 
independent of (r*/77*) in both quadratic and quartic models. 

Having focused on the end of inflation, previous analyses have neglected the contribution from the non-vanishing 
inflaton value to the gravitino mass. To estimate the significance of this effect, we also consider the evolution of the 
gravitino yield for the unrealistic case where = throughout inflation, and in Figure [Sl we show an example of 
our results for a quartic potential in both cases, with inflaton-independent gaugino masses. TJuring inflation the true 
yield is suppressed compared to the inflaton-independent gravitino yield due to the large gravitino mass supressing 
the collision term in Eq. ([3). We then observe a sudden increase in the true yield as the gravitino mass rapidly 
decreases and settles to its low-energy value, causing the m|/m^ term to dominate. The yield increases until the 
gravitino mass reaches rriQ^ and it is then just a matter of a few e-folds until the collision term becomes negligible 
and radiation fully dominates, at which point the yield freezes out. 




FIG. 3: The gravitino yield as a function of the number of e-folds for the quartic potential with inflaton-dependent (blue) and 
independent (red) gravitino mass, with inflaton-independent gaugino masses in both cases. In the left plot we vary and in 
the right plot we vary mg, 



In Figure |4j we show the difference in the thermal gravitino yield after freeze-out between inflaton-dependent and 
inflaton-independent gravitino masses. For large there is a negligible difference between the two cases. Increasing 
niQ^ reduces this difference, due to the rn^/tv?^ term never dominating, and in the large m^,^ limit Cq ^ . 




FIG. 4; Comparison between the thermally produced yield after freeze-out for inflaton-dependent (solid) and inflaton- 
independent (dashed) gravitino masses with a quartic potential and inflaton-independent gaugino masses. 

In Figure [Sj we plot the gravitino yield as a function of the number of e-folds, indicating where the above conditions 
are met. It is clear from this figure that the yield is not yet constant when pji = and that it changes quite 
drastically over a short number of e-folds until it becomes constant. Moreover, for the same parameters as in Figure 

t applying the standard reheating constraints using the temperature at which = or the temperature at which 
e yield freezes out results in ^ 10~^ and ~ 10"^*^, respectively, which are a few orders of magnitude lower 
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than the true yield. This is due to the cumulative effect of gravitino production throughout warm inflation and that 
has been neglected in earlier analyses of this problem. 




FIG. 5: The gravitino yield as a function of the number of e-folds for a quartic potential, indicating the number of e-folds at 

which slow-roll ends, pjj ~ p^, ~ and pji ~ exp(— 47Ve). These results correspond to = 2 x 10*, (Tt/i/«) = 1000 

and rrir, = 100 GeV. 
(jo 

Figure |6] shows the gravitino yield as a function of the number of e-folds for inflaton-dependent gaugino masses. 
We observe that, during inflation, gauginos are heavier than the gravitino and so the yield is larger than in the case 
where the gaugino masses do not depend on the inflaton fleld, once again due to the m|/m^ term in Eq. (j3|. The 
rise in Yq is due to mg/rriQ ^ {(f)/mp)^^ , so that as the inflaton fleld decreases this term enhances the yield until the 
gravitino mass settles at its low-energy value, w^^. As before, it is only a matter of a few e-folds until the collision 
term becomes negligible and the universe is in the radiation era. During this short number of e-folds the collision 
term evolves as Cq ~ and so the yield decreases until it freezes out. We also flnd that, as expected, the lower the 
value of JTigp , the larger the final yield is. 
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FIG. 6: The gravitino yield for inflaton-dependent gaugino masses as a function of the number of e-folds, for a quartic potential 
with Ctf, = 10* and different values of . 



C. Stable gravitinos 

Figure [t] shows the contribution of stable gravitinos to the current density parameter, fig/i^, for the quartic and 
quadratic potentials for inflaton-independent gaugino masses. Similarly, in Figure |8] we plot this contribution for 
inflaton-dependent gaugino masses. 

We can see that for sufficiently large C^, it is possible to avoid overclosure, VI. ^h? < 1, for a broad range of gravitino 
masses. This is related to the fact that increasing reduces the temperature of the radiation bath during inflation 
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FIG. 7: Contribution to the density parameter Q,Qh^ from an LSP gravitino for the quartic (left) and the quadratic (right) 
potentials, with inflaton-independent gaugino masses. Masses are given in GeV. 




FIG. 8: Contribution to the density parameter Q^h^ from an LSP gravitino for the quartic (left) and the quadratic (right) 
potentials, with inflaton-dependent gaugino masses. Masses are given in GeV. 

and hence reduces the thermal production. This can be achieved with lower values of in the quadratic model than 
the quartic, due to the former being flatter. If the gaugino masses depend on the inflaton, the overclosure problem 
becomes more severe. We can nevertheless satisfy the dark matter constraint for LSP gravitinos, fl^h^ < 0.1, if 
C4, ^ 1.5 X 10^" for the quartic and ^ 4 x 10^ for the quadratic potentials. At these large values of C^, there is 
little difference between inflaton-dependent and independent gaugino masses scenarios. 

For comparison with standard reheating predictions, we may define an effective reheat temperature as the tem- 
perature at which the gravitino yield becomes constant. In Figure [9] we illustrate the difference between the results 
predicted by Eq. ([t]) at this effective temperature with those obtained with the full numerical simulation for a quartic 
potential. 

We can conclude that, if the gaugino masses depend on the inflaton field, the standard reheating prediction is 
drastically different from the true warm inflation result, where the gravitino problem is more severe. If the gaugino 
masses are inflaton-independent then the standard prediction also leads to an underestimation of fl^h^. This implies 
that in warm inflation the effective reheat temperature needs to be somewhat lower than in standard reheating in 
order to avoid overclosure. For example, for a 1 TeV gravitino, standard constraints require Tr ^ 2 x 10^° GeV for 
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FIG. 9: Comparison between the predictions for Q,Qh^ using Eq. ([t]) at the efTective reheat temperature (dashed) and the full 
numerical simulation (solid) for inflaton-independent gaugino masses (left) and inflaton-dependent gaugino masses (right), with 
a quartic potential. All quantities are given in GeV. 



fig/i^ ^ 1, whereas in warm inflation we require Tr ^ 5 x 10^ GeV. 

As discussed earlier, in the cold inflation picture it is assumed that the yield of gravitinos at the reheat temperature 
is zero, Yq{Tii) = (see Section|ll|. This is perfectly valid in cold inflation, where due to the absence of a thermal bath 
during inflation, gravitinos are not produced. However, in warm inflation gravitinos are produced for the duration 
of inflation and so there is a non-negligible yield at the effective reheat temperature. Moreover, previous analyses of 
gravitino production during warm inflation assumed not only that the standard analysis at the end of inflation was 
applicable, but also that the gravitino yield froze out when — pn. We have seen that both these assumptions do not 
yield a good estimate for the gravitino abundance, both due to the cumulative effect of gravitino production during 
inflation and the fact that freeze-out does not occur until the universe if fully radiation-dominated, which occurs a 
few e-folds after inflaton-radiation equality. In particular, this results in an effective reheat temperature lower than 



previously estimated by more than one order of magnitude, as illustrated in Figure 10 
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FIG. 10: The temperature, in GeV, at which — pa and the effective reheat temperature, Tr, at which the gravitino yield 
freezes out, as a function of C4, for quartic potential. 

Constraints on the LSP gravitino also come from decays of the NLSP spoiling BBN predictions for light-element 
abundances. It is typical to assume that the NLSP is the MSSM-LSP and that it will only decay into the gravitino and 
Standard Model particles. The NLSP lifetime typically depends upon the gravitino mass and their mass difference. 
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"T-NLSP — 'TT'G; ^^'^ unless the gravitino is light and/or the NLSP is very heavy, it wiU be subject to BBN constraints. 
As the NLSP does not have planck-suppressed interactions with the other particles in the thermal bath, it may be 
in thermal equilibrium during inflation and freeze out in the radiation era. Its thermally produced yield will then be 
given by the freeze-out temperature, which places constraints on the NLSP and low scale gravitino masses but not 
on the warm inflation dynamics. In this respect the situation is the same as in cold inflation and, given that this is a 
model-dependent issue, we will not explore it any further, pointing the interested reader to the reviews in |20[ I21j. 



D. Unstable gravitino 



If the gravitino is the NLSP then we have the constraints from BBN on the primordial yield given in Section [IH 
and in order to obtain such low yields we must consider large values of C^. For « 100 GeV the bound on the 
gravitino-to-entropy yield is < 10~^^, which translates into ^ 10^^. Similarly, for m^, = ITeV the bounds are 
^ 10~^^ and Y? < 10~^^ for branching ratios into hadrons of Bh = 10^^ and Bh — I, respectively. This requires 

Ccf, ^ 10^^ and <^ 10^^, which are approximately the same for both the quartic and the quadratic potentials. 
Also, in the case of a gravitino NLSP, each_gravitino will then decay into one LSP. We can convert the primordial 

we show the hues at which ^l^sph^ 



gravitino yield into the LSP yield using Eq. fel). In Figure 
of values of in the quartic and quadratic potentials. 
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0.1 for a range 
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FIG. 11: Lines of flLSph'^ = 0.1 for various values of C,/, in the quartic (left) and quadratic (right) potentials, with NLSP 
gravitinos. The shaded region indicates where the LSP is heavier than the gravitino. Masses are given in GeV. 



We can see that the dark matter constraint can be satisfied for more reasonable values of Cd> than for the LSP 



gravitino. In particular, if tolsp 



100 GeV and m/^ 



1 TeV, the dark matter constraint is satisfied for Cs ^ 2.5 x 10^ 



(quartic) and ^ 6 x 10® (quadratic). If the NLSP gravitino mass, ^ 20 TeV, then it decays before BBN and 
the strongest constraint is given by the dark matter bound on the LSP. For m^^^ = 20 TeV the dark matter constraint 
is satisfied for ttilsp = 100 GcV with > 2.2 x 10^ (Q > 6 x 10®) and for mLSP = 1 TeV with > 2.5 x lO^'^ 
{Ccf, ^ 6 X 10^) in the quartic (quadratic) potential. 



CONCLUSION 



In this work, we have revisited the gravitino problem in warm inflation, focusing on thermal production which, 
providing the main difference from standard or cold inflation, places the strongest constraints on warm inflation 
dynamics. By performing a full numerical evolution of the gravitino yield into the radiation era we improve upon 
previous analyses. Firstly, in the context of thermal gravitino production, the effective reheat temperature is the 
temperature at which the gravitino yield freezes out and not the temperature at which the inflaton energy density 
equals the radiation energy density. This allows the temperature to drop by approximately an order of magnitude. 
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which lowers the final temperature at which gravitinos are produced compared to previous estimates. Secondly, we 
found that an analysis similar to standard reheating is in fact inadequate in describing gravitino production, due to 
the non- negligible yield produced throughout the whole duration of warm inflation. Finally, we have also taken into 
account the enhance particle masses during inflation due to supersymmetry breaking, in particular the gravitino and 
potentially the Standard Model gauginos. 

Taking all of these issues into account, our work shows, in particular, that the final gravitino yield is substantially 
lowered for stronger dissipative effects, as in practice this lowers the temperature of the radiation bath during warm 
inflation significantly. We have presented regions of parameter space where the LSP gravitino can satisfy the dark 
matter bound and, for an NLSP gravitino, we determined the regions where the LSP abundance does not exceed the 
amount of dark matter present in our universe and have given values of the dissipation parameter for which late 
decays do not spoil the predictions of BBN. 

Although thermal production is the dominant source of gravitinos during warm inflation, other non-thermal mech- 
anisms may play a role at a later stage. Gravitinos can, in particular, also be produced from particle decays, but 
due to the large Hubble parameter during warm inflation these decays will not take place until the radiation era, at 
which point the standard cosmological results can be used. They can also be produced from the direct decay of the 
inflaton field, although we have found that, in the sub-planckian regime, the dissipative ratio Q is necessarily large, 
which prevents the inflaton field from entering an oscillating phase. Figure 12 shows the inflaton field evolution as we 
artificially switch off dissipation at pr = p,f,, at which point oscillations immediately begin. 
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FIG. 12: Switching off dissipation at Pr = p^,, showing that the large dissipation keeps the inflaton fleld from oscillating in the 
radiation era. The dashed (solid) line corresponds to the case with (without) dissipation. 



Dissipation will actually switch off when the heavy fields are no longer kinematically allowed to decay into the 
light degrees of freedom, which depends on their low scale mass hierarchy. For example, if supersymmetry is indeed 
a solution to the gauge hierarchy problem, we may expect light scalar masses to lie close to the TeV scale and 
dissipation to switch off at temperatures of this order. Apart from these kinematical constraints, the form of the 
dissipation coefficient in Eq. (151 may actually hold down to very low temperatures. For example, to avoid exceeding 
the dark matter bound for the LSP gravitino, we require ~ 1.5 x 10^° and for 40 e-folds of infiation, if the coupling 
g 1, the system remains in the low-temperature regime down to T ~ lOMeV, at which point Pftjpcj, ~ 10^^. It is 
therefore unlikely in this case that any oscillations of the inflaton fleld may come to play a significant role in gravitino 
or, in fact, any entropy production. In particular, a significant dilution of the gravitino yield through a late inflaton 
decay along the lines proposed in [6j may be difficult to attain, although this may depend on the form of the inflaton 
potential, which goes beyond the scope of this work. 

Our analysis revealed that it is possible to satisfy the dark matter constraint for LSP gravitinos and LSPs produced 
from NLSP gravitinos at large values of the dissipation parameter C^, which requires large couplings and field 
multiplicites, pointing towards beyond the Standard Model scenarios. The gravitino problem is more severe for 
unstable gravitinos potentially spoiling the predictions of BBN, and in this case much larger values of are required. 

One should note that such large values of the dissipation coefficient are nevertheless required in order to overcome 
the severe eta-problem affecting monomial potentials for sub-planckian values. Above the Planck scale, the potential 
gets exponentially steeper with increasing field values, requiring larger values of to obtain 40-60 e-folds of inflation 



and also to suppress the resulting gravitino abundance. Using the full supergravity potential in Eq. ( 17) places a lower 
bound of > 10^ and ^ 2 x 10^ for 40 e-folds of infiation in the quartic and quadratic potentials, corresponding to 
(/>* ~ mp. Of course a non-canonical choice for the Kahler potential may alleviate this eta-problem, but supergravity 
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is in any case unlikely to be the complete theory near the Planck scale and so any analysis along these lines must 
be taken with a pinch of salt. It should nevertheless be emphasized that simple monomial potentials cannot yield 
the required number of e-folds for sub-planckian values without dissipation, which is an attractive feature of warm 
inflation despite the large field multiplicities and/or couplings required. 

One should bear in mind that observations may pose some constraints on the amount of dissipation present when 
the relevant CMB scales exit the horizon during inflation. In particular, an earlier analysis of non-gaussian effects on 
the primordial power spectrum showed that these depend logarithmically on the dissipative ratio at horizon crossing, 
Q*, placing a model-dependent upper bound on the parameter |63; . This analysis assumed, however, a constant 
dissipation coefficient, and more recently it was shown that for a generic T-dependence the non-gaussian parameter 
Jnl is largely independent of the value of Q* in the strong dissipative regime, yielding Jnl ^ C(10) within the 
observable window of Planck |53| . Hence, although the dynamics of second-order perturbations in warm inflation is 
not yet fully understood, we do not expect non-gaussianity to pose any significant constraints on our results. 

In this work, we have considered a general scenario where all the MSSM degrees of freedom are in thermal equilibrium 
during inflation. However, it has been pointed out in |60j that, in the low-temperature regime, mx ^ T, fermionic 
degrees of freedom may actually not thermalize, as both their contribution to the dissipation coefficient and their 
thermal scattering cross section are supressed compared to scalar fields. This is related to the structure of the 
superpotential ( 14 ) and the broken supersymmetry during inflation, which imply that the light fermions in the Y 
multiplets only interact via the heavy X bosons and fermions, whereas the light scalars have unsupressed interactions. 
Moreover, although the effects of gauge fields and their superpartners on the dissipation coefficient have yet to be 
analyzed in detail, their contributions to the dissipation coefficient may also be suppressed for sufficiently small gauge 
couplings. This would imply a thermal bath concurrent with inflation essentially composed of scalar particles, which 
would prevent gravitino production during inflation and eliminate the cumulative effect observed in our numerical 
simulations, at the same time requiring somewhat lower values of for sub-planckian inflation. Both fermionic and 
gauge degrees of freedom will nevertheless be 'reheaten' after inflation with either the exit from the low-temperature 
regime or the Hubble parameter dropping sufficiently in the radiation era. Although it requires further investigation, 
this may occur only at very low temperatures, as discussed above, in which case thermal gravitino production will be 
negligible. 

In cold inflation there is a tension between having a large enough reheat temperature for thermal baryogene- 
sis/leptogcncsis to occur and it being low enough to avoid overproduction of gravitinos and other unwanted relics 
(see e.g. [64 ). In warm inflation this can be alleviated, as a baryon asymmetry can be produced through dissipation 
itself [1]. Dissipation is an inherently out-of-equilibrium process, so the inclusion of baryon number and CP-violating 
interactions in the X and Y sectors in the superpotential Eq. ( 14 ) naturally leads to the production of a baryon asym- 
metry during inflation. In the low-temperature regime, the produced asymmetry is naturally small despite the large 
couplings and field multiplicities required for a sufficiently long period of accelerated expansion and, moreover, this 
may lead to distinctive baryon isocurvature perturbations in the CMB anisotropics spectrum that may be observable 
in the near future. Warm inflation thus exhibits several attractive features that address not only the problems of 
inflationary dynamics itself but also many of the associated cosmological puzzles. 

We would like to point out that the results from our analysis have a certain amount of crossover with cold inflation. 
Standard reheating is unlikely to be instantaneous and so the production of gravitinos will occur for the duration 
of the reheating phase. This will lead to an accumulated abundance similar to the one we have observed in warm 
inflation and so may change the standard reheating temperature constraints. The gravitino gets a mass from inflation 
and so, when the inflaton is oscillating about its minimum, the gravitino mass will also change at the same rate. If 
the oscillations are adiabatic, rh^j/mQ ^ Pscattcringj then this effect can be analysed for various potentials in a similar 
way to the analysis performed in this work. It may then result in significant differences in the thermal production of 
gravitinos during the standard reheating picture. 

With this work, we hope to have shed some light on gravitino production in warm inflation, with the way now 
paved for other potentials and dissipative coefficients to be analysed. In particular, the fact that inflation gives a mass 
to the gravitino may have a more signiflcant impact on the thermal production in other potentials. Our analysis also 
brought to light some issues that may be significant to standard reheating and we hope that this motivates further 
exploration of this topic. 
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